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Abstract 

We have studied the sensitivity to variations in the triple alpha and ^^C(a, 7)^^0 reaction rates, of 
the yield of the neutrino process isotopes ''Li, ^^B, ^^F, ^^^La, and ^^'^Ta in core collapse supernovae. 
Compared to solar abundances, less than 15% of ''Li, about 25-80% of ^^F, and about half of ^^^La is 
produced in these stars. Over a range of ±2(7 for each helium-burning rate, ^^B is overproduced and 
the yield varies by an amount larger than the variation caused by the effects of neutrino oscillations. 
The total ^^B yield, however, may eventually provide constraints on supernova neutrino spectra. 

PACS numbers: 26.30.Jk, 26.50. -hx, 14.60.Pq 
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About 10^^ neutrinos are emitted during a typical core collapse supernova explosion. 
For some time it has been known (see |l| for a detailed history) that interactions of these 
neutrinos with the stellar envelope can produce certain rare nuclei in abundances close to 
those observed in nature. These nuclei, called here the neutrino nuclei, include ''Li, ^^B, ^^F, 
i38La, and ^^OTa [l|,3. 

It was pointed out ^ that the production of some of the ^^"^Ta and most of the ^^^La by 
the neutrino process was sensitive to the electron neutrino temperatures, and might serve to 
probe the value of the neutrino oscillation parameter sin^2^i3. Recently, js-S] showed that 
the yields of ^Li and ^^B in supernova explosions are also sensitive to sin^ 29is and to whether 
the neutrino mass hierarchy is normal or inverted. In both cases, this sensitivity arises 
because neutrino oscillations can change the neutrino spectra produced during core collapse 
supernovae, increasing the average energies of the and i>e and affecting the synthesis of the 
neutrino nuclei. Since two of the main goals of neutrino physics 6| are to determine better 
the value of sin^ 2^13 and the nature of the mass hierarchy, the possibility that the observed 
abundances of the neutrino nuclei might constrain these quantities is of great interest. 

Their use for this purpose depends, however, on the robustness of the stellar yield pre- 
dictions. Studies of the dependence of nucleosynthesis on the helium burning reaction rates 
have shown 0, H] that both the yields of the more abundant nuclides and stellar structure 
are significantly affected. Since the neutrino nuclei result from neutrino induced spallation of 
abundant progenitor nuclei, their production depends on the abundances of these nuclei and 
on their location within the star, and thereby on the rates of the helium burning reactions. 

In this paper, we examine the changes in the production of ^Li, ^^B, ^^F, ^'^^La, and ^^°Ta 
caused by changes in the astrophysical helium burning rates within their uncertainty limits, 
and compare the yield changes of ^Li, and ^^B, with the predicted js-S] effects of oscillations. 

We then discuss how, and whether, the neutrino process nuclei can be used to constrain 
the neutrino spectra from supernovae. We find that the constraints provided by neutrino 
process nucleosynthesis are interesting but not yet definitive. Because of the great interest 
in these issues it appears that a major effort to sharpen these constraints is warranted; a 
discussion of important measurements and calculations is given below. 

We used the KEPLER code j9-12] to model the evolution of 15, 20, and 25 solar mass 
stars from central hydrogen burning up to core-collapse; a piston placed at the base of the 
oxygen shell was then used to simulate the explosion. Following ^ we assumed a total 
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energy of 5 x 10^^ ergs per neutrino species, i.e. a total of 3 x 10^^ ergs energy release in 
the supernova explosion. Mass loss processes were included. The neutrino spectra were 
approximated by Fermi-Dirac distributions with a zero degeneracy parameter, a luminosity 
exponentially decaying after onset of core collapse with a time-scale of 3 s and a constant 
neutrino temperature: T = 4 MeV for Uf, and Ue', T = 6 MeV for u^, u^, v^-^and Vr- For 



further details, see 



0.aa 



13|. These choices are consistent with estimates of neutrino 
emission intensity and time dependence from supernovae fl^ . 



Initial stellar abundances were taken from both Anders & Grevesse 



151 ] and from Lodders 



16|, hereafter AG89 and L03. The L03 abundances for C, N, O, Ne are roughly 15%- 



25% lower than those of AG89, whereas the abundances of heavier elements are roughly 
!-5% higher. For calculations of neutrino process cross sections we used the results from 
Briefly, in that paper the charged and neutral current cross sections were first used to 
calculate the excitation spectra of the product nuclei; experimental data and O^w shell-model 
estimates were used to determine the Gamow- Teller response for ^^C and ^°Ne (leading to 
^^B and ^^F) and RPA estimates for the J < 4 multipoles for all other transitions. The 
SMOKER statistical model code was then used there to follow the ensuing decays. The 
7 process contributions to the yields are also included in our calculations, but are important 
only for ^^°Ta. 

Isotope yields were stored at nine key points of stellar evolution js], Q. As anticipated, 
^Li and ^^B were produced essentially only during the supernova stage; their yields are 
shown in Fig. 1. Here an initial "A" (or "L") label means that the calculations were done 
for the AG89 (or L03) abundances, and a final "A" (or "C") means that the ^2C(a,7)^^0 
(or triple-alpha) rates were varied by ±2cr from their central values. These central values 
were, resp., 1.2 times the rate recommended by Buchmann [l^ with cr = 25% and that 
recommended by Caughlan and Fowler 171 with o = 12%. For the ^^C(q!,7)^^0 rate, the 

p n 

central value is that commonly used in calculations with the KEPLER code [3, 119[; it is 
consistent with recent measurements [2^. The energy dependence obtained by Buchmann 
was used for all calculations. The labels also give the stellar mass, or 3 STARS, the average 
for the 15, 20, and 25 Mq stars using a Scalo 2l| Initial Mass Function (IMF) with a slope 
of 7 = —2.65. For normalization purposes we compare to the production factor for oxygen. 
Since ^^O is made mainly in massive stars, a production factor ratio near one is consistent 
with all (or most) of an isotope being made in a primary neutrino process (as are ^Li, ^^B, 
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FIG. 1: Production factors of ^Li and ^^B compared to those of ^^O for various reaction rates. 
The left hand column shows the results when the triple-alpha reaction rate R^a is varied about 
the central value of 1.0, the rate of ref. The right-hand column shows the results when Ra,i2 
is varied about the central value of 1.2. The value for ''Li has been multiplied by a factor of 10. 
An example of the range of variation in ^^B yield predicted in 31] is shown as a band in the upper 
right-hand panel. The dotted line at 0.4 is the production factor ratio that would give the solar 
abundance of ^^B not made in the galactic cosmic rays. For more information see the text. 

and i^F). 

Examining first the results for average production in the three-star sample, and assuming 
that this is a reasonable approximation of the total production process, we see that only 
10%-15% of ^Li is made in the neutrino process. This is not surprising, since there are 



many processes, including the Big Bang, that make or destroy ^Li and that are not fully 
understood. On the other hand, for most values of the reaction rates ^^B is overproduced, 
even if one ignores production by cosmic rays. 

Fig. 2 shows the results for ^^^La. and ^'^'^Ta. Here we show also results for the pre-SN 
stage, the time when the contraction speed in the iron core reaches 1000 km sec~^, since 
production during that stage is not negligible, especially for ^^°Ta. A detailed examination, 
however, shows that most of the ^'^^La and ^^°Ta that is ejected in the SN is not what was 
present in the pre-SN stage; that is mostly destroyed by the SN shock and most of what is 
ejected was newly synthesized during the explosion {2!. Ill|. 

Since these two isotopes are secondary products (produced from pre-existing spallation 
targets) a production ratio of about two for ^^^La and ^^°Ta (see the dotted line on Fig. 2) 
would be necessary to reproduce the solar abundance. An additional complication is that 
our models do not distinguish production in the short lived ground state from that in 
the long-lived isomeric 9~ state in ^'^^Ta; a better, but still approximate, treatment 22 1 
gives an isomer production of about 40% of the total production. It then appears that the 
production of ^^'^Ta is roughly consistent with the solar abundance, given the uncertainties in 
the production calculations, and that the production of ^^^La corresponds to about half the 
solar abundance. ^^F is a primary product, and it appears that 25%-75% of solar ^^F could 
be made by the neutrino process. This complicates the determination of the importance of 
other sources such as AGB stars and Wolf-Rayet winds. 

We now consider whether a comparison of the observed abundances of ^Li and ^^B to SN 
model predictions can place constraints on the neutrino oscillation process, as was suggested 
in Refs. 3|-[5|. These investigations were for a 16.2 Mq star, using parameters almost identical 
to those we have used, except that the explosion energy, T^^, and Tp^ were 1.0 Bethe, 3.2 
MeV, and 5.0 MeV instead of 1.2 Bethe, 4.0 MeV and 4.0 MeV. The near equality of the 
average neutrino energies should yield similar production for the two models in the absence 
of neutrino oscillations. 

In js-S] neutrino oscillations produce significant increases in ^Li production, up to 75%, 
as sin^ 26*13 increases from 10^^ to 10^^ for the normal neutrino hierarchy-the changes are 
much smaller, around 15% for an inverted hierarchy. The changes for ^^B are small for 
either hierarchy, around 20%. The number ratio A^(''Li )/A^(^^B) is assumed to be less 
susceptible then the absolute yields, to systematic uncertainties in the calculations and has 
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FIG. 2: (Top left): Production factors (three star averages) for ^^'^La and ^^''Ta using Anders 
and Grevesse abundances and varying the triple alpha rate. (Top right): Same but varying Ra,i2- 
(Lower left): Same for 15 Mq star, varying Ra,i2- (Lower right): Production factors (3 star 
average for^^F. The solid curves are for the final abundances and the dashed curves are for the 
pre-SN stage. All production factors are ratios to that of ^^O. 

approximately 50% changes for the normal hierarchy, and less than 10% for the inverted 
hierarchy. This provides, in principle, some hope that observation of an enhanced ratio 
could place a lower limit on the value of sin^ 2^13 and eliminate the option of an inverted 
hierarchy. 

The variations with the helium burning reaction rates, of the production of ^Li by the 
neutrino process are relatively small, about 20%, but it will be difficult to untangle the 
relatively small amount of neutrino-produced ^Li from other sources of ^Li. One might hope 
that observations of ^Li in pre-solar grains would make it possible to isolate the effects of 
individual supernovae. Unfortunately, few, if any, relevant observations have been made to 
date. Lithium isotopic ratios have been measured in very large SiC grains by Gyngard, et 
al. 23|, but Li is very volatile and is not expected to condense into SN SiC grains 24 1. 
Thus, while its production does not depend strongly on the helium burning rates, other 
considerations limit its usefulness. 



The situation for ^^B is also unclear. Hoppe et al. |25| measured B in supernova SiC 



grains. Their measured isotopic ratio, ^^B/^°B = 3.46 ± 1.36, is, however, consistent with 
laboratory contamination by solar system B (^^B/^'^B= 4.045). They conclude that at most 
30% of the measured B is attributable to the neutrino process. Their Fig. 5 25| shows 
a corrected value reduced by a corresponding factor of three. After this correction, the 
abundance of B is lower than expected, by over an order of magnitude. They consider some 
possible explanations, but it appears that the grain formation process is not well understood. 

Moreover, as we see in the top four panels of Fig. 1, variations in ^^B yields with reaction 
rate are large, a factor of two or more, making the uncertainties in its predicted ratio to 
^Li and in its absolute value much larger than the effects predicted by 3|-l5|. We conclude 
that this approach to constraining neutrino oscillations will not be productive until there 
are significant improvements in the helium burning rates discussed here, as well as in grain 
observations and their interpretation. 

Before considering whether the gross production rates might eventually provide a con- 
straint on neutrino spectra, we need to examine the total production of ^^B. Both ^^B and 
^"^B are made in the galactic cosmic rays (GCR) with the ratio ^^B/^*^B lying between 2.2 



and 2.5 



26 



271 ]. The ±0.15 uncertainties in the ratio reflect, main 



cosmic ray sources and the propagation model. (For reviews see 



y, uncertainties in the 
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observed meteoritic ratio 4.045 



29[). We take as the 



16l ]: the other recent abundance summaries 



15 



30| quote 

results within 0.5% of this value. Since we find that the neutrino process makes little ^"B 
^iig^iog gg^^ about 42 ±4% of the ^^B must be made in the neutrino process. A possible 
contribution of (so far unobserved) low energy cosmic rays would increase the ^^B/^'^B ratio 
in cosmic ray production 26| but it has been found 3l| that this process is energy inefficient 
and unlikely to produce a significant amount of ^^B. The neutrino process ^^B should then 
be compared to about 0.4x (solar ^^B)-this comparison is made in Fig. 1. 

Summarizing, except for ^Li and ^^F, which have other known production sites, it appears 
that production by the neutrino process (and partially, for ^®°Ta by the gamma process), 
as shown in the 3-Star panels of Figs. 1 and 2, is within a factor of three of the observed 
abundances except for extreme values of the rates. 

It has been shown previously that neutrino process yields increase strongly for larger 
neutrino energies. For example, increasing the temperature of v^, and neutrinos from 4 



to 6 MeV increases the yields of the neutrino nuclei by factors from 1.5 to 2 



2j . Similar 



changes were obtained for z/^. This strong dependence raises the possibility of constraining 
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the ranges of allowable neutrino temperatures, spectral shapes and neutrino intensity. Such 
constraints depend on the robustness of the model predictions, and thereby on the nuclear 
rates, on the neutrino interaction cross sections, on the form of the neutrino spectra, and 
on the astrophysical modeling uncertainties of the underlying stellar models. It is probable 
that the best limits will be obtained for ^^B. The neutrino interaction cross sections for 
^^C can be more reliably calculated than those for the heavier nuclei, because the strong 
Gamow- Teller cross sections are mainly experimentally based, and shell model estimates can 



replace RPA calculations for the L > cross sections |32|-|34| . The constraint imposed by 
the meteoritic and GCR ^^B/^'^B ratios is also useful in determining the appropriate SN ^^B 
yield. 

Taken at face value, it seems that significantly harder neutrino spectra than we have 
used are improbable-the yield is already overestimated. But improvements in the neutrino 
process calculations are necessary to make this constraint credible. It appears likely 35 1 
that the uncertainty in the triple alpha rate will be halved in the near future and there are 
major efforts to improve the ^'^Cj a, ^Y ^O rate. Better estimates of neutrino spectra can also 
be employed; it is now known 36|, |37| that the mean energies of the various neutrino species 
are more similar than had been thought, and that the high energy tail is suppressed by 
inclusion of inelastic scattering processes. The mean energies and second energy moments 
of these new spectra are, however, similar to those of the Fermi-Dirac distributions we have 
used, differing by less than 12 % in all cases-the second moments are related to the neutrino 
process cross sections. The astrophysical model uncertainties are a remaining difficulty, but 
these should also be reduced using techniques informed by 3-D calculations. 

To summarize, we explored changes in the core-collapse supernova yields of ^Li, ^^B, 
i9p^ ^^^La, and ^^°Ta that arise from changes in the triple alpha and ^^C(a,7)^®0 reaction 
rates within their ±2cr uncertainties. We found that the rate changes result in factor of two 
changes in the production of ^^B in a 15 Mq star. This, for the present at least, rules out 
the techniques proposed 3|-l5| to constrain the neutrino oscillation parameter sin^ 2^13. For 
the assumed neutrino spectra there is significant overproduction of ^^B for all values of the 
rates we have used (±2(t). It seems reasonable to expect that a factor of two improvement 
in the precision of neutrino-process nucleosynthesis can be achieved, especially for ^^B. This 
may provide a constraint on the neutrino energy spectrum. If one assumes that model 
calculations can accurately fix the spectral shape, neutrino process nucleosynthesis could 



8 



provide an estimate of the neutrino flux from supernovae and a check on supernova models 
that does not depend on occurrence of (infrequent) supernova explosions. 
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